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Abstract 



Context. X-ray emitting features originating from the interaction of supernova shock waves with small interstellar gas clouds are 
revealed in many X-ray observations of evolved supernova remnants (e.g. Cygnus Loop and Vela), but their interpretation is not 
straightforward. 

Aims. We develop a self-consistent method for the analysis and interpretation of shock-cloud interactions in middle-aged supernova 
remnants, which can provide the key parameters of the system and the role of relevant physical effects like the thermal conduction, 
without the need to run ad-hoc numerical simulations and to bother of morphology details. 

Methods. We explore all the possible values of the shock speed and cloud density contrast relevant to middle-aged SNRs with a set 
of hydrodynamic simulations of shock-cloud interaction, including the effects of thermal conduction and radiative cooling. From the 
simulations, we synthesize spatially and spectrally resolved focal-plane data as they would be collected with XMM-Newton/EPIC, an 
X-ray instrument commonly used in these studies. 

Results. We devise and tune up two diagnostic tools, the first based on the mean-photon energy vs. count rate scatter plot and the 
second on the spectral analysis of the interaction region, that can be used to highlight the effects of thermal conduction and to derive 
the shock speed in case of efficient conduction at work. These tools can be used to ascertain information from X-ray observations, 
without the need to develop detailed and ad-hoc numerical models for the interpretation of the data. 
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1. Introduction 

Supernova remnants (SNRs) are known to be a privileged lab- 
oratory to investigate the physical and chemical evolution of 
the galactic interstellar medium (ISM) and the mass distribution 
of the plasma in the Galaxy. Multi -wavelength observations of 
evolved SNRs (e.g., iGraham et 10 ^995: Bocchino et al. 2000; 



iPatnaude et al.l2002l: 



Nichols & Slavin,.2004;,Miceli et al...2005i) 



can be a useful tool to investigate the physics of SNRs, for 
instance the interaction of the remnants with inhomogeneities 
(clouds) of the ISM. However, this interaction involves many 
non-linear physical processes (e.g. radiative losses and thermal 
conduction) which make the analysis of the observations quite 
difficult. A further limitation comes from the superposition of 
5_i different emitting regions along the line-of-sight (hereafter LoS) 
Ci ^ and, in most cases, the data interpretation is not unique. 

A powerful approach in the data analysis is based on hy- 
drodynamic and MHD simulations of the shock-cloud inter- 
action, which takes into account the most relevant physical 
mechanisms (e.g. thermal conduction, radiative cooling, etc.), 
and on the comp arison of th e model results w ith observations . 
Previous studies dOrlando e t al. 2005, 2006; Miceh et al.ll2006t 
lOrlando et al.l 2008*) were devoted to investigate, through numer- 
ical modeling, the interaction of SNR shock fronts with small in- 
terstellar gas clouds. The scope includes: i) to investigate the role 
of the different physical processes at work on the dynamics and 
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energetic of the shocked cloud, and ii) to analyze accurately the 
SNRs observations through their comparison with model results. 

As a part of this project, we already investigated the role 
of thermal conduction and radiative cooling on the evolution of 
the shocked cloud in the unmagnetized limit. We explored two 
physical regimes i n which each of the two physical processes 
in turn dominates (lOrlando et al.l l2005l hereafter Paper I) and 
found that, in general, the thermal conduction determines the 
evaporation of a fraction of the shocked cloud, forming a hot 
and tenuous gas phase (the corona) surrounding the cloud core. 
In the presence of an organized interstellar magnetic field, the 
thermal conduction is known to be inhibited across the mag- 
netic field lines and the radiative cooling can be enhanced due to 
magnetic plasma confinement. We explored the role played by 
the magnetic-field-oriented thermal conduction and the radiative 
cooling during the shock-cloud interaction, considering differ- 
ent configurations of the magnetic field (Orlando et al. 200^. 
We found that the magnetized cases fall in between the limit 
of completely suppressed thermal conduction and the unmagne- 
tized limit with conduction. 

Our numerical models have also been used to make predic- 
tions on the expected X-ray emission from the shock-cloud in- 
teraction. We showed that the X-ray emitting structures do not 
trace the morphology of the flow structures originating from 
the shock-cloud interaction and that the shocked clouds are 
visible more easily during the early phases of their evolution 
(Orlando et al. 2006, hereafter Paper II). 

However, the big effort done in the modeling of shock-cloud 
interaction and its X-ray emission has not been counterbalanced 
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by a rigorous methodology in the comparison between X-ray 
observations of SNR shells and models. The high resolution in- 
struments on board XMM-Newton and Chandra have provided 
us with excellent images and spectra of SNRs which are always 
much more complicated than the ideal cases treated in numerical 
simulations. Therefore, a straightforward comparison between 
models and observations is still difficult, and this tends to ham- 
per our understanding of the details of the physical processes 
which are at the ba se of the X-ray radiation from SNR shells. 
iMiceli et al.l (12006 ') made one of the first attempt in filling the 
gap between models and observations. They compared the X- 
ray observations of an isol ated knot in the northern rim of the 
Vela SNR (Vela FilD; Mic eli et al.ll2005h with an ad-hoc hydro- 
dynamic model; the comparison showed that the bulk of the X- 
ray emission in the knot originates in the cloud material heated 
by the transmitted shock front, but significant X-ray emission 
is also associated to the cloud material which evaporates in the 
intercloud medium, under the effect of the thermal conduction. 
While this strategy has proved to be winning, it is quite model- 
dependent, in the sense that it is based on an accurate and strict 
morphological and spectral comparison which may be time and 
resource consuming. 

On the contrary, the idea behind this paper is that the explo- 
ration of the parameter space of the shock-cloud model already 
performed in Paper I and Paper II, along with the extension pre- 
sented here toward still unexplored values of the cloud density 
contrast, may be used to devise a quick and effective method- 
ology for the interpretation of current generation X-ray satel- 
lite observations of shock-cloud interactions, without the need 
of running ad-hoc numerical models. Our intent is to provide 
easy-of-use recipes that allow to extract from the data many of 
the key parameters governing the evolution of shocked clouds, 
by comparison with a set of model quantities normalized in a 
way to eliminate the dependence from unnecessary details, like 
e.g. the exact morphology of the hit cloud. In particular, our 
scope includes to devise a diagnostic tool able to quickly as- 
sess if the spectral results obtained in the interaction regions are 
dominated by thermal conduction, a physical effect whose con- 
tribution to the X-ray emission is modulated by the magnetic 
field and, therefore, still in debate and uncertain. 

The paper is organized as follows: Sect. |2] briefly describes 
the numerical setup, the physical parameters of the problem, and 
the method to synthesize, from the numerical simulations. X-ray 
observations as they would be obtained with X-ray observato- 
ries; Sect. [3] presents the results of the numerical simulations; 
in Sect. |4] we describe the diagnostic tools devised in this pa- 
per and apply the methods, as an example, to X-ray observations 
reported in the literature; in Sect.|5]we draw our conclusions. 



2. Hydrodynamic modeling 

We model the three-dimensional interaction of a SNR shock 
front with an ISM cloud in the same way we have done in Paper 
I, to which the reader is referred to have more details. We sum- 
marize here the main model features. The cloud is assumed to 
be small compared to the curvature radius of the shoclfl and in 
pressure equilibrium with the unperturbed isothermal and homo- 
geneous ambient medium; we consider, therefore, a planar shock 
front and an isobaric cloud, spherical for simplicity. The shock 
propagates with a Mach number Al » 1 in the ambient medium. 
The post-shock initial conditions are given by the strong shock 



Table 1. Parameters of the simulated shock-cloud interactions. 
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" Shock Mach number. * Density contrast cloud/ambient medium. 

Velocity of the SNR shock. Temperature of th e post-shock ambient 
medium. '' Cloud crushing time dKlein et aljl994h . 



limit jZel'dovich & Raizeij 19661) . The fluid is assumed fully ion- 
ized, and is regarded as a perfect gas. 

The plasma dynamics is described by solving numerically 
the time-dependent fluid equations of mass, momentum, and en- 
ergy conservation (see Eqs. 1-5 in Paper I). The m odel takes into 
account the thermal conduction (ISpitzeii 119621) and the radia- 
tive losses from an o pticall y thin plasma (e.g. Ravrnond & SmithI 
|T977[ iMewe et all [T985i and iKaastra & Mewd \266(h . The 



' This assumption is valid for a 1 pc cloud in the middle-aged SNRs 
Vela and Cygnus Loop, whose shell has a radius > 10 pc. 



thermal conduction includes the free-streaming limit (satura- 
tion) on flie hea t flux (Cowie & McKee 1977, Giuliani 1981 
Bor kowski et al.l [1989, Fadevev et alj l2002t and references 
therein). Our calculations also include a passive tracer associated 
with the cloud material to trace its motion during the evolution. 
A discussion of the assumptions of the model and their influence 
on the results is presented in Sect. 14.41 

The numerical code is flash (iFrvxell et al . 1 120001) . a multi- 
dimensional hydrodynamics code for simulating astrop hysical 
plasmas, which uses the paramesh (iMacNeice et al.l2000l) library 
for block-structured adaptive mesh refinement (AMR), and has 
been customized with numerical modules that treat thermal con- 
duction and optically thin radiative losses (see Paper I for de- 
tails). The initial configuration, the boundary conditions, and the 
AMR setup of the simulations used here are the same as those 
adopted and discussed in Paper I. 

We consider, as a reference case, the M - 50 shock model 
described in Paper I; we then explore the parameter space by 
varying, alternatively, either the Mach number, M, or the density 
contrast cloud/surrounding medium, x- In the reference model 
(RCmSOclO), the unperturbed ambient medium is at tempera- 
ture Tism = 10^ K and particle number density Wism = 0.1 cm"-', 
the spherical isobaric cloud has a radius = 1 pc and den- 
sity contrast =10 (particle number density Wci - A'"ism - 1 
cm"^). The SNR shock front is planar at Mach number Al = 50 
and temperature Tp^h = 4.7 MK. In the other simulations, the 
Mach number varies in the range 40 < Al < 60 (correspond- 
ing to shock temperatures in the range 3 MK < Tp^h < 7 MK) 
and the cloud density contrast in the range 3 < < 30 (cor- 
responding to particle number density of the cloud in the range 
0.3 cm"-' < Hci < 3 cm"^). We note that in Paper II, we have 
already partially explored the variation induced by a different 
choice of the shock speed (we have considered Al = 30 and 50 
atx - 10). Here, we present for the first time the results for dif- 
ferent density contrasts. These ranges are representative of most 
of the shock-cloud interaction regions observed in evolved SNRs 
(e.g. Vela, Cygnus Loop, and G296.5-H10.0). 
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Figure 1. x ~ ^ parameter space (adapted form Paper I). The 
lines are derived for length-scale / = 1 pc, and for an unper- 
turbed ambient medium with temperature Tism = 10^ K and par- 
ticle number density Wism = 0.1 cm"^ (see Paper I for details): 
the solid line separates regions dominated by radiative cooling 
(on the left) from regions dominated by thermal conduction (on 
the right); the dashed line marks the density contrast, ^, above 
which the cooling time-scale T,ad is shorter than the cloud crush- 
ing time Tec; the dotted line marks the values x above which 
the thermal conduction time-scale Tcond is shorter than Tcc- The 
parameter pairs explored are marked with crosses. 

The effects of thermal conduction on the shocked cloud evo- 
lution have been fully investigated in Papers I and II. Since we 
are interested in deriving some diagnostic to be used in real X- 
ray observations, we have compared runs with this physical pro- 
cess (hereafter RC runs) with other runs without it (hereafter 
HY runs). As shown in a previous work (Orlando et al. 2008), 
shock-cloud interactions in an organized interstellar magnetic 
field fall in between these two limits (i.e. HY and RC cases). 
A summary of all the simulations discussed in this paper is in 
Table [U while Fig. [T] shows the simulations in the;^' - M. param- 
eter space. As discussed in Paper I (cfr. Fig. 2 in Paper I), this 
plot can be used to evaluate if radiative cooling is competitive 
with respect to thermal conduction for a given run. For example, 
the shock transmitted into the cloud is strongly radiative in runs 
RCm40clO and RCm50c30. On the other hand, the thermal con- 
duction dominates over the radiative losses in all the other cases 
(i.e. RCm50c03, RCmSOclO and RCm60clO): the cloud is ex- 
pected to evaporate on a time-scale comparable (in RCm50c03 
and RCmSOclO) or shorter (in RCm60clO) than Tcc- 

2. 7. Synthesis of the X-ray observations 

The output of the numerical simulations is the evolution of 
temperature, density, and velocity of the plasma in the spa- 
tial domain. From the density and temperature values, we syn- 
thesize spatially and spectrally resolved X-ray observations 
with the XMM-New ton/EPIC-pn X-ray imaging spectrometers 
( IStruderetalj2001h . The method can be easily extended to other 
X-ray instruments. 

The emission measure in the yth domain cell is calculated as 
emj = «Hj^J' where n^^ is the particle number density in the 
cell, and Vj is the cell volume. We assume that the direction 
of the LoS corresponds to the y axis (in the Cartesian coordi- 
nate system), perpendicular to the direction of propagation of 



the SNR shock front, and that the depth along the LoS is 10 pc 
(a typical value for the shell of evolved SNRs, such as in Vela 
and in Cygnus Loop). We then derive the distributions of emis- 
sion measure versus temperature, EM(r), integrated along the 
LoS for each {x, z), in the temperature range 5 < log r(K) < 7 
(divided into 50 bins, all equal on a logarithmic scale). From 
the EM(r) distributions, we synthesize maps of X-ray emission 
and X-ray spectra, using th e MEKAL spectral synthesis code 
dMewe et all 1 19851; iKaastrai 1992, Kaastra & Mewe 2000), as- 
suming solar metal abundances (Grevesse & An ders.. 1991 ). 

We assume the source to be at a distance Dsnr = 500 pc 
(as, for instance, in the case of Cygnus Loop) and we filter 
the spectra thro ugh an ISM absorption column density, A^h = 
5 X 10 cm omson &McCammonl 119831) . according to 
typical values derived fro m SNR observations at that distance 
(e.g. Patnaude et al.ll20d2l) . The absorbed X-ray spectra are then 
folded through the instrumental response to obtain focal plane 
spectra. The exposure time is assumed to be t ^xo = 10 ks 
for EPIC-pn (see, for instance, iMiceh et al.ll2005|) . The photon 
counts are randomized in each energy instrumental channel of 
the focal-plane spectra according to Poisson statistics, using the 
rejection method applied to the Poisson distribution (iPress et aP 
.1986) . X-ray emission maps are produced in selected energy 
bands, assuming a spatial resolution of 4 arcsec; the X-ray im- 
ages are convolved with the corresponding point spread function 
(PSF), as given by Ghizzardi (2002) for EPIC-pn. 

The final products are X-ray simulated observations, spa- 
tially and spectrally resolved, in a format virtually identical to 
that of real observations collected with EPIC-pn. To such data, 
we apply the standard methods of analysis commonly used for 
X-ray observations. 

3. Results 

3.1. Light-curves 

The detectability of the shock-cloud collision in the X-ray band 
is expected to depend on the M and x parameters. From our 
simulations, we derive the X-ray light curves of the region as- 
sociated to the shocked cloud, to understand at which stage of 
the interaction the visibility of the cloud is maximum. Such a 
region is selected in each synthesized EPIC-pn count rate im- 
age in the [0.3 - 2.0] keV band (typically selected for the anal- 
ysis of evolved SN R shock-cloud interaction; see, for instance, 
IMiceh et al.l2005h . by considering all the pixels having a median 
energy of X-ray photons, & dHong et alj |2004| ). which is less 

than 90% of the & derived for the surrounding mediurrQ; from 
these pixels, then, we evaluate the average counts s"' per pixel. 
Fx, normalized to the value derived for the intercloud medium. 
The X-ray light curves reported in Fig.|2]show when the shocked 
cloud is detectable (when the normalized Fx is higher than 1) 
and its X-ray luminosity is maximum. 

In general, we find that the higher is M, the higher the nor- 
malized Fx at each stage of the evolution (see Fig. |2]l. In most 
cases, the shocked cloud is visible in the X-ray band in the in- 
terval 0.1 < f < 1.5 Tcc. The thermal conduction makes the 
shocked cloud brighter than in cases without conduction, broad- 
ening the peak in the X-ray light-curve for any Mach number. In 
fact, the conduction contributes to the cloud heating, increasing 
the amount of cloud material above 1 MK and emitting in the X- 
ray band. Note also that the conduction makes the shocked cloud 

- In fact, the shocked cloud material is expected to be cooler than the 
surrounding medium. 
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Figure 2. X-ray light-curves of the shock-cloud interaction re- 
gion in the [0.3 - 2] keV band (see text for details on the def- 
inition). The time dependent surface brightness is normalized 
to the average value of the post-shock intercloud region. Upper 
panel: runs without thermal conduction (HY runs) with;^' =10 
and M = 40, 50, 60. Middle panel: runs with thermal conduc- 
tion (RC runs) with;^' =10 and M = 40,50,60. Lower panel: 
RC runs with M = 50 and;^- = 10, 30. 



hardly detectable in cases with x ^ ^ due to the quick evapo- 
ration of the cloud. For instance, in run RCm50c03, there are 
no pixels with S below 90% of the value derived for the back- 
ground and no light-curve can be reported in Fig.|2] On the other 
hand, the conduction makes shocked clouds with ;^ > 20 (whose 
evolution is strongly dominated by the radiative losses) partially 
visible during the very early phases of the evolution (see lower 
panel in Fig.|2]l; these clouds would be not detectable in X-rays, 
in the absence of thermal conduction, being their temperature 
Tgci ~ 2.5 X 10^ K. In these cases, only the thermally conducting 
corona is detected, being the core much cooler than 1 MK. 



3.2. Spectral analysis 

The shocked cloud is detectable with EPIC-pn during the early 
phases of the shock-cloud interaction (f < 1.5 Tcc; see Fig. |2]l 
as a bright knot surrounded by a diffuse region (see right pan- 
els in Fig. O. We focus the spectral analysis in the interval 
0.4 Tcc ^ t < 1.4 Tcc when the shocked cloud is the bright- 
est in all the models. For each sampled X-ray image, we select 
spatial sub-regions in the computational domain and analyze the 
X-ray spectra extracted from each of them. To select spectrally 
homogeneous regions, we derive maps of the median energy of 
X-ray photons, £, from the EPIC-pn data that allows to con- 
vey at the same ti me both sp atial and spectral information on the 
emitting plasma jHong et al. 2004). Since the shocked cloud is 
cooler than the surrounding medium, the median photon energy, 
S , of the bright region is lower than that of the surroundings (see 
left panels in Fig. O. Thus we select subregions with a median 
photon energy 0.5 keV < S < 0.6 keV (to identify the knot), 
and subregions with 0.6 keV < £ < 0.7 keV (for the diffuse 
region, DR). The knot corresponds to the brightest portion of 
the X-ray image in Fig. [3] and the DR selects the intermediate 
brightness region surrounding the knot (compare left and right 
panels in Fig. [3]). 

It is important to stress that this definition of the extraction 
regions for spectral analysis are completely independent from 
the morphology of the X-ray emission. This has the great ad- 
vantage that it can be straightforwardly applied to any current 
X-ray telescope observation for which the mean photon energy 
map can be computed. Moreover, it makes the spectral analysis 
independent, in first approximation, from the details of the shape 
of the ISM clouds, which in reality may be much more complex 
than the ideal spherical cloud proposed in our model. 

The extracted spectra have a total number of photons rang- 
ing between 10"* and 4 x 10^, adequate for a detailed spec- 
tral analysis. The focal plane spectra hav e been analyzed us- 
ing the spectral fitting package XSPEC (Arnaud 1996) and ap- 
plying a multi-temperature fit to each spectrum. All the ex- 
tracted spectra are well fitted with two MEKAL components of 
an optically-thin thermal plasma in collisional ionization equi- 
librium (Meweetal. 1985, Kaastra & Mewe 2000), with so- 
lar abundances, and filtered thro ugh the interstellar absorption 

orrison & McCammonlll983l) . We have applied this proce- 
dure to any model listed in Table [1] and in the interval 0.4 Tcc < 
f < 1.4 Tcc in step of 6t - 0.1 Tcc. We present here, as an exam- 
ple, the results obtained in the reference models HYm50clO and 
RCm50clO. 

Figure |4] shows the temperature, T, and the emission mea- 
sure per unit area, em = EM/Areg (where A,eg is the area of 
the selected region), of the isothermal components fitting the 
EPIC-pn spectra. When the thermal conduction is completely 
suppressed (run HYm50clO), the spectra of both the knot and 
the DR at the different epochs are described, in general, by two 
isothermal components with temperatures Tio„ ~ 1 MK and 
Thigh ~ 4.5 MK; the emission measure of the hot component is 
emhigh ~ 5 X 10"^ cm"^ in all the spectra, whereas emiow ranges 
between 3 x lO'*^ and 8 x lO'*^ cm . Note that the temperature 
of the hot component is close to the temperature of the shocked 
ambient plasma Tpsh ~ 4.7 MK, whereas Tiow is slightly higher 
than the temperature of the shock transmitted into the cloud, 
r,ci ~ 0.8 MK (see Paper I). 

In run RCm50clO, the spectra are again described by two 
isothermal components, but with some differences due to the 
thermal conduction. In particular, Tiow is higher and T^gh is 
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Figures. Median energy maps (left) and EPIC-pn count rate 
images (right) in the [0.3 - 2.0] keV band derived for run 
RCmSOclO at the three labeled times during the evolution. The 
pixel size is ~ 4 arcsec and the exposure time is 10 ks. The im- 
ages are smoothed with a boxcar of width cr - 12 arcsec. The 
white contours mark the cross-section of the cloud on the plane 
of the image, identified by zones consisting of the original cloud 
material by more than 90%; the contours superimposed to the 
median energy maps mark the bright knot (red) and the diffuse 
region (diffuse region, DR; black). 



lower than the values derived in HYmSOc 1 0, being the difference 
larger for the DR than for the knot. Also emio„ is systematically 
lower and emhigh higher than the values derived in HYmSOc 10. 
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Figure 4. Best-fit values of temperature (upper panel) and emis- 
sion measure per unit area (lower panel) for the EPIC-pn spec- 
tra extracted from the bright knot and from the diffuse region 
(DR) in runs HYmSOc 10 (blue) and RCmSOclO (red) at differ- 
ent epochs between 0.4 and 1 .4 t^c- The dashed lines in the upper 
panel mark the temperatures expected for the shock transmitted 
into the cloud (Tsd ~ 0.8 MK) and for the shocked ambient 
plasma (Fp^h ~ 4.7 MK). 



In general, we find that the cold component is the most sensi- 
tive to the thermal conduction, showing the largest differences 
between HYmSOc 10 and RCmSOclO. The origin of these differ- 
ences is explained by comparing the results of the spectral fitting 
with the distributions of emission measure per unit area versus 
temperature, em(r), from which the extracted spectra originate. 

Figure |5] shows these em(r) distributions together with the 
results of the spectral fitting. In general, the distributions of both 
the knot and the DR are bi-modal. The cold peak around T » 2 
MK is due to the shocked cloud gas; the hot peak at T ^ S MK is 
due to the shocked ambient plasma surrounding the cloud. The 
best-fit values are localized around these maxima and, therefore, 
can be associated to the shocked cloud gas (cold component) and 
to the shocked ambient plasma (hot component). 

In HYmSOclO, the em(r) distributions of both the knot and 
the DR do not change significantly during the evolution (except 
at f ~ 1.2 Tec, when the shocks transmitted from the front and 
from the rear of the cloud interact; see Paper II), with the two 
bumps steadily centered around the temperatures expected for 
the shock transmitted into the cloud (Tsd ~ 0.8 MK) and for 
the shocked ambient plasma {Tpsh ~ 4.7 MK), respectively. By 
comparing HYmSOclO with RCmSOclO, the main effects of the 
thermal conduction are: i) to smooth the em(r) distributions, be- 
cause of a transition region formed between the inner part of the 
cloud and the ambient medium in which the density decreases 
and the temperature increases smoothly in the radial direction 
(see Paper II); and ii) to shift the first bump to higher temper- 
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Figures. Distributions of emission measure 
per unit area versus temperature, em(r), in the 
range of temperature 5 < log r( K) < 7 for the 
knot (solid) and for the diffuse region (DR, dot- 
ted) at the labeled times during the evolution. 
The figure shows the em(r) derived from run 
HYmSOclO (on the left) and from RCmSOclO 
(on the right). The figure also shows the results 
of the spectral fitting to the EPIC-pn spectra ex- 
tracted from the selected regions (symbols as in 
Fig-Ell- The errors are at 90% confidence level. 



atures due to the gradual thermalization of the shocked cloud 
material to the temperature of the shocked ambient plasma (see 
also Paper II). 

As a result of the conduction effects, the amount of plasma 
above 1 MK increases in RCmSOclO, making the shocked cloud 
brighter in the X-ray band (see Fig.|2|. Also, the changes in the 
em(r) distributions due to the thermal conduction determine the 
differences in the results of the spectral fitting for HYmSOclO 
and RCmSOclO (see Fig. |4|; for instance, the shift of the first 
bump in em(r) to higher temperatures leads to higher Tiow, and 
the smoothing of em(r) leads to lower emiow in RCmSOclO. 
Note also that the effects of the thermal conduction are the 
largest in the em(r) distribution of the DR, being the plasma 
of the corona surrounding the cloud core subject to efficient heat 
conduction. As a consequence, the cold fitting component in the 
DR is, in general, hotter than that in the knot (see also Fig. |4|, 
whereas the opposite is true in HYmSOclO (i.e. the temperature 
of the shocked cloud material is never higher than the tempera- 
ture of the shock transmitted into the cloud). 

4. Diagnostics 

4.1. Median energy vs. count-rate scatter plot 

Since the thermal conduction modifies the temperature and den- 
sity structure of the shocked cloud (see Paper I), its effects may 
be expected in the comparison of & maps (related to the spatial 
distribution of temperature) with count rate maps (related to the 



spatial distribution of mass density). We deri ve, therefore, S ver - 
sus count rate scatter plots (see, for instance. iMiceli et"al]|2005h : 
we first divide the range of count rate [0.01 - 0.20] cnts s"' into 
100 bins (all equal on linear scale); then, from the EPIC-pn count 
rate images in the [0.3 - 2.0] keV band, we derive the median 
photon energy of all the pixels belonging to the same count rate 
bin. Fig. |6] shows the scatter plots derived for HY (upper panel) 
and RC (middle and lower panels) runs at selected epochs be- 
tween 0.4 Tec < f < 1 .4 Tec (when the shocked clouds are visible; 
see Fig.|2]l. All these plots are characterized by a clear descend- 
ing trend and, then, in most of the cases, by a much flatter fall 
(cold plateau): the higher the count rate, the lower is the me- 
dian energy and, therefore, the lower is the average temperature 
along the LoS. The descending branch and the cold plateau are 
the signature of the shock-cloud collision: the former roughly 
coiTesponds to the DR and the latter to the bright and cold knot 
defined in Sect. 13.21 

In HY runs, scatter plots derived for;^' - 10 and different Ai 
show a similar shape, characterized by a very steep descending 
branch and a well defined cold plateau (see upper panel in Fig. 
|6|l. The descending branch shows an abrupt transition between 
the intercloud material (highest £ ) and the shocked cloud ma- 
terial (lowest fi). In RC runs (see middle panel in Fig. |6ll, the 
thermal conduction makes the slope of the descending branch 
flatter than that of HY runs, being the absolute value of the 
slope smaller for higher Ai. The flattening of the descending 
branch reflects a smooth temperature and density structure of the 
shocked cloud due to the heat conduction that leads to the grad- 
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Figure 6. Median photon energy, &, versus count rate scatter 
plot derived from the EPIC-pn data in the period 0.4 Tcc < f < 
1 .4 Tec- Upper panel: runs without thermal conduction (HY runs) 
with;if = 10 and M - 40,50,60. Middle panel: runs with ther- 
mal conduction (RC runs) with - 10 and M - 40,50,60. 
Lower panel: RC runs with = 50 and ;y = 3, 10,30. The 
dashed lines in the upper and middle panels show the slope of 
the £ , versus count rate scatter plot derived from the analysis of 
EPIC data of Vela FilD (iMiceli et al.ll2005h . 



ual growth of a transition region from the inner part of the cloud 
to the ambient medium (see Paper I). Note that the scatter plots 
of runs with different and same Al = 50 are virtually indistin- 
guishable (see lower panel in Fig.|6ll independently of the role of 
radiation or conduction. This is due to the fact that the contribu- 
tion to X-ray emission invariably comes from regions dominated 
by thermal conduction. In particular, in run RCm50c30, the X- 
ray emission originates from the thermally conducting corona, 
being the cloud core at temperatures T < 2 x 10^ K. 
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Figure?. The figure compares the temperature and the emis- 
sion measure per unit area of the cold component derived for 
the knot with those derived for the diffuse region (DR). The blue 
stars (red crosses) mark the values derived in run HYm50clO 
(RCm50clO) and reported in Fig.|4l the black symbol with error 
bars mark the values derived from the analysis of EPIC data of 
Vela F ilD (regions 4, knot, and 1, diffuse region, in lMiceli et alJ 
I2OOI . 



In summary, these scatter plots can be very useful to deter- 
mine the role of the thermal conduction through the slope of 
the descending branch and to derive hints about the speed of 
the shock (in case of efficient conduction at work). On the other 
hand, the plots are poorly useful to infer the actual density con- 
trast of the cloud. 



4.2. Temperature and emission measure ratios 



The spectral analysis discussed in Sect. l3.2l suggests that the cold 
fitting components describing the knot and the DR are sensitive 
to the thermal conduction; we propose, therefore, to use them 
as a diagnostic tool to trace the efficiency of conduction. Fig. [T] 
compares the temperature and emission measure values derived 
for the knot with those derived for the DR (and reported in Fig.|4| 
in runs HYm50clO and RCm50clO. The runs with/without ther- 
mal conduction are clearly separated in the plot, the thermal con- 
ductive case being localized in the bottom-right quadrant and 
the pure hydrodynamic case in the top-left quadrant. This result 
is determined by the development of the thermally conducting 
corona in RCm50clO and, therefore, is expected to be general. 
In particular, as discussed in Sect. 13.21 the thermal conduction 
smooths the first bump in the em(r) distributions and shifts it 
to higher temperatures, being this effect larger for the DR than 
for the knot. As a result, the cold fitting component derived 
for the DR is, in general, hotter than that derived for the knot 
(i.e. 7'iow[DR] > riow[knot]), whereas the opposite is true (i.e. 
riow[DR] < riow[knot]) if the conduction is suppressed. Also, 
the smoothing of the em(r) distribution due to the conduction 
is the largest for the DR (because the corona surrounding the 
cloud core is subject to efficient heat conduction), leading to the 
smallest values of emiow (see also lower panel in Fig.lUl. 

In summary, the temperature and emission measure ratios are 
an excellent way to determine the role of the thermal conduction 
in the evolution of the system. 
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4.3. An example of model vs. observation comparison 

Our study shows that evidence of thermal conduction at work 
during the shock-cloud interaction may be found in the spectral 
analysis of X-ray data. As discussed in the previous sections, 
the & vs. count-rate scatter plot and the temperature and emis- 
sion measure ratios can be efficient diagnostic tools to derive 
the shock speed and the role of thermal conduction, which, in 
turn, is linked to the magnetic field configuration, as shown by 
lOrlando et al.l (l2008l) . In this section, we challenge the above di- 
agnostic tools and show that they can be easily used in the anal- 
ysis of X-ray data, by comparing our model results with X-ray 
observations reported in the literature. 

In particular, we focus on a well-studied region, the "FilD 
region", that is an isolated, bright X-ray knot in the northern 
rim of the Vela SNR. Because of i ts proximity (~ 250 pc, 
iBocchino etaP 1 19991 ICha et al.l 11999 ) Vela is an ideal target 
for this kind of study, allowing us to observe the interaction 
of the SNR shock fro nt with relatively sm all clouds, like FilD 
(~ 2 X 10''^ cm; see iMiceU et all l2005h in great det ail. The 
analy sis of an XMM-Newton observation of FilD (Mice li et al.l 
l2005h has shown that its X-ray spectra can be modeled by an 
optically-thin plasma with two thermal components (at ~ 1 MK 
and ~ 3 MK, respectively) with inhomogeneous volume distri- 
butions along the line of sight. The cold component dominates in 
the brightest region that is surrounded by a diffuse region with 
harder X-ray emission. To interpret these results, iMiceli et al.l 
(|2006|) developed a detailed hydrodynamic model of FilD, syn- 
thesized X-ray emission maps and spectra from the model, and 
compared them with the data. Their analysis has shown that 
the X-ray and optical emission of FilD can be explained as the 
result of the interaction of a SNR shock (with Mach number 
Ai = 57) with an ellipsoidal cloud 30 times denser than the 
intercloud medium ; the estimated interaction time is ~ 0.32 Tcc. 
iMiceli et al.l 12006) proved that the two components originate in 
the cloud material heated by the transmitted shock front and by 
heat conduction between the cloud and the hotter, shocked in- 
tercloud medium. FilD, therefore, is an ideal benchmark for our 
model, being a case for which the thermal conduction has been 
proved to be at work. Since the parameters used in our simula- 
tions are slightly different from the parameters deduced from the 
observations (including the shape of the cloud), we do not ex- 
pect a perfect match, but the comparison will nonetheless give 
us many useful information. 

Among the runs presented here, the one matching the den- 
sity contrast of the shock-cloud interaction is RCm50c30 (see 
Tab. [TJ. As akeady discussed, the shocked cloud with such a 
density contrast (x ~ 30) would not be detectable in X-rays if 
the thermal conduction is suppressed, being its estimated tem- 
perature Tsci ~ 3.5 X 10^ K. On the other hand, inspecting Fig.|2] 
we note that, in run RCm50c30 (dashed line in the lower panel), 
the shocked cloud is visible in X-rays for a short time inter- 
val (0.1 - 0.6 Tcc) and has the largest surface brightness at the 
evolutionary stage estimated for FilD (~ 0.32 Tcc). We expect a 
brighter emission for higher values of Ai (see middle panel in 
Fig. |2|i. As shown by our simulations, the detected X-ray emis- 
sion originates in the cloud material dominated by thermal con- 
duction, confirming the relevance of conduction in the evolution 
of FilD. 

IMiceli et al.l (l2005l) analyzed the spectra extracted from the 
knot and the DR composing the FilD region. Thus, we can derive 
the temperature and emission measure ratios of the cold com- 
ponents derived by these authors and plot them in Fig. Q The 
observed values lie in the bottom-right quadrant of the figure. 



confirming once again that in FilD the thermal conduction is ef- 
ficient, in perfect agreement with the independent conclusion of 
Miceli et al. (2006). The diagnostics in Fig. Qis of easy imple- 
mentation and we suggest it as a standard to check for the role 

of conduction. _ 

IMiceli et al.l (l2005l) derived also a £ versus count rate scat- 
ter plot for FilD that can be compared directly with the corre- 
sponding scatter plots derived with our models. We overplotted 
a best-fit power-law model (with ind ex = -0.25) to the d ata of 
the FilD region reported in Fig. 4 of IMiceli et al.l ( l2005b . con- 
sidering the count rate as a free parameter, a reasonable choice 
if we consider that the actual value of the count rate depends on 
the actual LoS extension, which is poorly known. The slope of 
the observed scatter plot is rather flat and cannot be reproduced 
by models without the thermal conduction (see upper panel in 
Fig. O. On the other hand, the observed slope is reproduced 
quite well by our RC runs, in agreement with the evidence that 
the thermal conduction plays an important role in the evolution 
of FilD. Al = 60 seems to be the model which best reproduces 
the slope, in very good agreement with the value obtained by 
iMiceli et al.. (.2005.) (M = 57) with a detailed analysis. 

4.4. Limits of the model 

In our simulations, we parametrize the thermal conductivity us- 
ing the classical Spitzer's conductivity and the saturation limit, 
assuming essentially laminar thermal conduction in the all spa- 
tial domain. However, regions of strong turbulence of differ- 
ent strength and extent can develop in the system (especially 
in shock-cloud interactions dominated by radiative cooling), for 
instance at the shear layers along the cloud boundary or at the 
vortex sheets in the cloud wake. The turbulence in these regions 
may have a significant effect on the thermal conduction, leading 
to significant deviations of thermal conductivity f rom i t s laminar 
values (e.g. lKulsrudll983l:lNaravan & Medvedevl200ll ; lLazarianl 
2006). As a result, the thermal conduction may be inhomoge- 
neous due to the presence of turbulence. On the other hand, the 
deviations of thermal conductivity from its laminar values are 
expected to be relevant in the shocked intercloud medium, thus 
not affecting our main conclusions on the effects of thermal con- 
duction on the shocked cloud and on the applicability of the di- 
agnostics devised here. 

In our model, we do not account for the possible effect of 
the back-reaction of accelerated cosmic rays on shock dynam- 
ics. In the case of high Mach number shocks, a part of the shock 
power may be dissipated into cosmic rays acceleration, resulting 
in the increase of the shock compression ratio. The distribution 
function of non-thermal particles and the bulk flow profile in 
the shock upstream region are sensitive to the total compression 
ratio. Thus, even a moderate efficiency of particle acceleration 
may reduce the p ost-shock ion and electron temperatures (see, 
e.g. Eq. 18 in Bvk ov et"an i2008). with implications on the X- 
ray emission. This effect is expected to be large for shocks with 
high Mach number (as, for instance, in young SNRs), but not 
in middle-aged SNRs (to which this paper is focussed on) for 
which no non-thermal emission has been detected. Even if par- 
ticle acceleration were not negligible, the effect relevant to our 
diagnostics would only be to slightly reduce the efficiency of 
the thermal conduction because of the low er post-sho ck temper- 
ature. As in the case of magnetized clouds |Orlandoet al. 2008j), 
shocked clouds with considerable particle acceleration would 
fall in between the limit of completely suppressed thermal con- 
duction (HY runs) and the unmagnetized limit with conduction 
(RC runs) discussed in this paper. 
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Recently [Pittard et aP (l2009h have demonstrated that the tur- 
bulence plays an important role in shockcloud interactions, and 
that environmental turbulence adds a new dimension to the pa- 
rameter space. In particular these authors have shown that the 
turbulence is mainly generated around the cloud boundary and 
in the cloud wake after ~ t^c', the main effect is that clouds sub- 
ject to a highly turbulent post-shock environment are destroyed 
significantly quicker than those within a smooth flow: the larger 
the cloud density contrast the higher is the effect of turbulence 
(for instance, for ~ 100, the effect of the post-shock turbulence 
dominates in the shock-cloud interaction). On the other hand, an 
efficient thermal conduction makes the cloud boundary smooth 
very quickly (see Paper I), and the turbulenc e grows more slowl y 
around clouds with a smooth density profile dPittard et al.ll2009l) . 
Thus we are confident that our results are valid and the diag- 
nostics proposed here can be applied to the clouds with moder- 
ate density contrast (i.e. the effects of turbulence poorly influ- 
ence the shock-cloud interaction at early evolutionary stages for 
t < Tec) and smooth density profiles (i.e. the growth of turbu- 
lence around clouds is slow), that we considered here. 

Finally, our model does not account for the incomplete 
electron-ion temperature equilibration in the post-shock region. 
Equilibrium may not be complete early during the shock-cloud 
interaction (t 0.1 Tcc). In that phase emission models includ- 
ing non-equilibrium should be applied and the initial part of the 
light-curves presented in Fig.|2]may change. On the other hand, 
the diagnostics discussed in this section refer to the shock-cloud 
interaction at 0.4 Tcc < f < 1.4 Tqc, when the hypothesis of tem- 
perature equilibration ca n be considered realis tic for the shock 
velocities explored here (iRakowski et al.ll2003l) . 



simulations. While this method cannot be considered a substitu- 
tion for a detailed approach to the study of shock-cloud interac- 
tions based on the developing of ad-hoc (and time-consuming) 
numerical models, nevertheless, we have demonstrated that the 
diagnostic diagrams we have presented can be very useful to de- 
rive some of the parameters of the system and the role of ther- 
mal conduction in a very quick and straightforward way. These 
information can be used in turn for a more detailed model, if 
necessary. 

The method can be applied to imaging X-ray observa- 
tions of middle-aged thermal SNR shells (like Vela or Cygnus 
Loop), as for instance those obtained by the XMM-Newton 
and Chandra X-ray satellites. We used as a benchmark the 
XMM-Newt o n/EPI C observations of the Vela Fil D region of 
'Mic eli et alj (l2005h . from which, independently, Miceli et"al] 
(200^^ have found strong evidence of thermal conduction at 
work during the shock-cloud interaction, using a detailed ad-hoc 
numerical model. We found that our method is quite effective in 
recovering quickly the shock speed and the effects of the thermal 
conduction. 
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5. Summary and conclusion 

In a series of previous paper (Paper I, Il. lOrlando et al.l200 8'). we 
have investigated the X-ray emission arising from the interaction 
of SNR shock waves with isolated gas clouds with the scope of 
identifying the plasma structures that mainly contribute to X- 
ray emission detectable with current X-ray instruments. In this 
work, we extend our parameter space considering clouds with 
different value of density contrast and devise diagnostics in the 
X-ray band revealing the shock speed, which is one of the fun- 
damental parameter governing the shock-cloud interactions, and 
the cloud evaporation under the effect of the thermal conduction. 

In particular, by performing a series of spectral fittings to the 
simulated data of the shock-cloud interaction region, we proved 
that there are at least two interesting diagnostic diagrams that 
can be used: 

- the median energy vs. count rate scatter plot (Fig.|6l), which 
gives a direct estimate of the shock speed and a hints about 
the effects of the thermal conduction; 

- the temperature and emission measure ratios between the 
knot and diffuse region of the cloud (Fig. |7]l, which gives 
a direct estimate of the role of the thermal conduction in the 
evolution of the system. 

We stress that the regions (the knot and the DR) which must 
be selected to derive the diagnostic diagrams are defined entirely 
on the basis of mean photon energy maps, and not on the ac- 
tual shape of the X-ray emission. Therefore, the method is very 
well-posed and independent, in first approximation, to the ac- 
tual shape of the ISM clouds, which could be more complex 
that the ideal spherical cases considered in our hydrodynamic 
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